The extraplanar stellar populations of highly-inclined disk galaxies IC 2233, IC 5052, 
NGC 4631 and NGC 5023 are analyzed with the goal to quantify their vertical extent and 
structure. Based on the single-star photometry, we separate different stellar populations and an- 
alyze their spatial distribution. On archival images obtained with the Hubble Space Telescope 
ACS/WFC the surroundings of these galaxies are well resolved into stars with the red giant 
population (RGB) identified far above the galaxy mid-planes. We find that there are a profound 
change in slope of the number density profile of the evolved RGB stars at extraplanar height 
of 4 - 8 kpc, that possibly reflects a truncation of the thick disk and reaching the 2-3 times 
more extended oblate stellar component ("the halo"). This structure is consistent with our pre- 
■ ' vious studies of both edge-on and face-on disk galaxies in the Local Universe (Tikhonov et al. 

' 2005a,b) and allow us to improve the spatial model of the stellar components of a typical spi- 
CsJ , ral galaxy. In NGC 4631, the revealed asymmetry of its stellar thick disk and halo, is likely 
, caused by the neighbor dwarf galaxy NGC 4627. Based on the tip of the red giant branch method 
^ ' (TRGB) we estimated a distance of 10.42 + 0.38 Mpc for IC 2233, 5.62 + 0.20 Mpc for IC 5052, 
7.11 + 0.13 Mpc for NGC 4631 (6.70 + 0.15 for its satellite NGC 4627) and 6.14 + 0.15 Mpc for 
[ NGC 5023. We confirm the presence of slight extraplanar metalUcity gradient of evolved stars at 



stars. 



<^ NGC 5023 and IC 5052, based on the systematic changes of the colour distribution of red giant 

m 
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Abstract. 

1. Introduction 

Understanding the structure of the stellar component of disk galaxies is an important step for un- 
raveling their formation and evolution. Although challenging to measure, the extraplanar stellar 
content around galaxies is a potentially-important probe of the most ancient stellar populations, 
as well as of the role of that different external and internal processes play in the evolution of 





Abadi et alJ2003 ; Brook et al J2OO4 ; 


Bullock & Johnsonl2005 ): disk heatine (e.s... 


KrouDa 


2002 


GnediJ2003 


; or the ejection of of gas from the underl^ 


fing disk, which may cool so that 
). 


stars can form (e.g.. 


Howk & Savaeell999 


: Savaee et alJ2003 



The stellar outskirts around our own Galaxy and some other large spiral galaxies has been 



extensively s 
them (e.g.. 


tudied over the last few years and revealed faint stellar envelopes in many of 


SacketJ 


1999 


( 


Iluillandre et al. 


1999; 


van der Marell 


2002; 


Fereuson et al. 


2002; 


Dalcanton & Bernsteinll2002 


Zibetti & Fereusonl2004 


). The most extensive and detailed study 



of the outer stellar populations of the Milky Way and M 31 revealed complex patchy struc- 
tures with numerous streams, which appear to be the remnants of disrupted companion galaxies 



* Based on observations with the NASA/ESA Hubble Space Telescope, obtained at the Space Telescope 
Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., 
under NASA contract NAS 5-26555. 
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(e.g., 



Brown et alJl2003 



: IZucker et al.l2004l) . However, the chemical patterns of the stars of the 
dwarf spheroidal (dSph), and irregular (LMC and SMC) satellite galaxies are predominantly dis- 
tinct from the stars in each of the kinematic components of the Galaxy. This result rules out con- 
tinuous merging of low-mass ga laxies similar to the currently observed satellites during the for- 



Smecker-Hane et alJll998 



McWilliam & Smecker-Hanel2005 



mation of the M ilky Way (e.g 

Venn et al]l2004 ). At the same time, the detection of extraplanar cold and ionized gas regions at 



large distances from the mid-planes of some spiral galaxies, which must be powered by hot stars 
suggests that while most of the outer, elusive stellar components are predomina ntly evolved, 
some fraction of these stars are not a single burst ancient stellar popul ation (e.g., 
Howk & Savagell999 



mi 



HoonesetalJl999t 



iTiillmann et 31120031) 



Keppel et al. 



In the literature, the definition "halo" and "thick disk" is rather unclear, as well as the method 
of their separation are not well established yet. While some of the recent studies of the outer 
regions of spiral galaxies rely on multicolour surface photometry, successful detection is only 
possible in the deep ob servation s of the edge-on galaxies with a prominent thick disk and hal o 



Dalcanton & Bernsteini2002t iPohlen et al. 



2004 



Zibetti et alJl2004 



Zibetti & Fergusonl2004 ). 



For example. 



Dalcanton & Bernstein 



(I2002h analyzed the sample of 47 edge-on spiral galaxies 
and have found a presence of the thick disk in 90% of them, based on the results of multi- 
colo ur surface photomet ry. By stacking of 1054 rescaled images of edge-on galaxies in the Sloan 



DSS 



Zibetti et al 



i2004l) detected an extended low surface brightness halo, however the "mean" 
colours of this component (g - r = 0.65 and r - i - 0.6) cannot be easily explained by normal 
stellar populations, and may suggest contamination by ionized gas emission. The important other 
result obtained at surface photometry studies is that both thin and thick disks are truncated, with 
mean scaleheight and scalelength of the thick disk is several times larger than that of the thin 
disk (Pohlen et al. 2004). Giv en that there hav e also been some no n-detections of luminous thick 

r lprvet alj[l999 : 



disk or halo components (e.g[ 



Zheng et al 



1999ft . the question arises about the 



ubiquity, origin and further evolution of these outer structures. 

The alternative technique of tracking out resolved stars to explore various galactic subsystems 
sidesteps the difficulties inherent in quantifying extremely faint diffuse emission. The measure- 



stellar components in Milky Way and in a 


handful of nearest galaxies ( 


Gilmore & Reid 


1983 


Chiba & Beer ^200(1 


Prochaska etalJ200n 


Worthev et alJ2005tlIbata et alJ2005l). For the more 



distant galaxies this information is unavailable, and a set of the model dependent assumptions 
need to be considered. Aside from the spatial distribution law and kinematics, stars from each 



taining considerable age and chemical information i \ 


/allenari et alJboOO 


Prochaska et alJ2000 


Chiba & Beers 


2000; 


Saraiedini & Van Duvne 


2001 


WilliamJ 


2002; 


Brewer & CarnevI 


2004) 



On the basis of multi-colour single-star photometry, one is able to separate different stellar pop- 
ulations, and suggesting that each of the components has a different dominant stellar population, 
to constrain their spatial distribution. 



4 Tikhonov et al.: Stellar outskirts around edge-on galaxies 

Table 1. The ACS/WFC Survey Fields 



Galaxy Region Date Band Exposure ID N, 



IC2233 


SI 


2004-04-29 


F814w 


2x350 


9765 


27120 






2004-04-29 


F606w 


2x338 


9765 




NGC4631 


SI 


2003-08-03 


F814w 


2x350 


9765 


115315 






2003-08-03 


F606w 


2x338 


9765 






S2 


2004-06-09 


F814w 


2x350 


9765 


108980 






2004-06-09 


F606w 


2x338 


9765 




IC5052 


SI 


2003-12-13 


F814w 


2x350 


9765 


63108 






2003-12-13 


F606w 


2x338 


9765 




NGC5023 


SI 


2004-07-02 


F814w 


2x350 


9765 


49628 






2004-07-02 


F606w 


2x338 


9765 





Inspired by the unparalleled angular resolution and sensitivity of the HST WFPC2 and ACS, 
we began a project to analyze the outermost stellar structure of spiral galaxies. In related works, 
we have analyzed the the spatial distribution and properties of the outer stars in the nearby face- 
on (M 81 and NGC 300), as well as edge-on (NGC 55, NGC 891, NGC 4144, and NGC 4244) 



spiral galaxies jTikhonov et al.ll2005 



Tikhonov & Galazutdinoval2005l) . In the most massive of 



the studied disk galaxies we have found a profound change in slope of the number density pro- 
file of the evolved red giant branch (RGB) stars, that likely reflects a truncation of the the thick 
disk and reaching the pure halo stellar component. The similar change of the number density 



and M 82 iTikhonovl2005l) . These results strongly support the ubiquitous presence of the trun- 



122^ 



gradient of RGB stars were also reported for lower mass spiral M SSjj^uillgnyr^^I 
Rowe et al.l2005 ), as well a s in two massive irregular galaxies: IC 10 ( Drozdovskv et al J2003 ) 



cated thick disks in the majority of disk galaxies, as well as large extend of the stellar outskirts 
in all directions. 

In this paper, we selected four high-inclination nearby disk galaxies for an in-depth study of 
their extraplanar steUar structure: IC 2233, IC 5052, NGC 4631 and NGC 5023. Our primary 
objectives are to: (i) quantify the stellar population variations associated with extreme outer sub- 
structure, and (ii) derive stringent constraints on the age and metallicity of stars in the far outer 
disk. The choice of high-inclination galaxies ensures minimal contamination by stars belonging 
to the dominant galactic disk component, and provides the more robust separation of the thick 
disk and halo stellar populations. 
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Fig. 1. The DSS-2 images of the galaxies with the ACS/WFC footprints superposed. The edge of 
thick disk and halo extent are shown by the ellipses. 

2. Observations and stellar photometry 

To study the resolved stellar population of galaxies, we obtained ACS/WFC images available 
in the HST archive. Digital Sky Survey images of these galaxies with ACS/WFC footprints 
are shown in Fig. [l] An interior and exterior ellipses delineate an estimated edges of the thick 
disk and halo. The observational data are listed in Table where ID is the HST program num- 
ber and ^stars is the number of detected stars . Images were reproces sed through the standard 



ACS STScI pipeline, CALACS, as described by 



Pavlovskv et al 



1 20051) . All data have first order 



bias subtraction, dark subtraction, and bad pixel masking applied. Images are then corrected for 
the flat-field response. We combined images using a latest version of the Multidrizzle pack- 
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age ( iKoekemoer et alJ2002h . which provides an automated method for distortion-correcting and 
combining dithered images. Multidrizzle also corrects for gain and bias offsets between WFC 
chips, identifies and removes cosmic rays and cosmetic defects. We generate point-source cata- 
logs fr om the stacked direct image using an automatic star-finding program DAOPHOT/FIND in 
MIDAS I Stetsonll994 ), applying a 3. Oct per pixel detection threshold. The measurements of their 
magnitudes were conducted via point-spread-function (PSF) fitting that is constructed from the 
isolated 'PSF-stars' as it is implemented in DAOPHOT/ALLSTAR. The background galaxies, unre- 
solved blends and stars contaminated by cosmetic CCD blemishes were elimi nated from the final 



fists, using their peculiar ALLSTAR characteristics, \SHARP\ > 0.3, x > 1-3 (IStetsoiIl994 . We 
derived an aperture correction from the 2.5 to the 10 pixels (C'S) aperture radius by analyzing 
profile of the PSF-stars. 

The F6Q6W and 14W instrumental magnitudes have been transformed to standard magni- 
tudes in the Kron-Cousins system adopting our empirically-derived transformation: 



(y-/)= 1.321 -(v- + 1.133 
/ = / + 0.059 ■(¥-!) + 25.972 



(1) 



These transformations have been calculated using the photometry of 60 stars with < (V - 
/) < 5.0, detect ed in both, ACSAVFC and WFPC2 images of irregular galaxy IC 10, and the 



prescriptions of 



Holtzman et al 



lll995h for the WFPC2 data. The photometric uncertainties of 
this calibration are dominated by the accuracy of the PSF-photometry of these stars with a median 
value of 0'."03 in /-band and 0?'04 in V. 
In 



Tikhonov et al 



1 20051) and 



Tikhonov & Galazutdinoval 1 120051) we applied these equations 



to the single-star photomefi-y of the ACS/WFC images of NGC 55 and NGC 4244. Based on 
the derived magnitude of the tip of the red giant branch method, we estimated distances to these 
galaxies, which agree well with measurements obtained with WFPC2 and other instruments. 



With the extinction law of lCardelli et al 



tinction, adopting the foreground extinction of 



1 1989 ) we correct all the photometry for Galactic ex- 



Schleeel et al 



1 19981) . and assuming that internal 



extinction in all galaxies is negligible in their outer ("halo") parts. 



3. Results 

3.1. CM diagrams 

The results of photometry for all stars detected in the galaxies are presented as colour-magnitude 
diagrams (CMD) in the Fig. |2 while Fig. |3] show the CMDs of the stars from the outermost 
extraplanar areas The features of these diagrams resemble those of spiral galaxies. All the 
CMDs are characterized by a noticeable red plume, while strength of the blue plume varies a lot. 
The spatial variations are apparent from the varying strengths of the blue and red plumes, and 

' Note, that outer part of NGC 4631 also contain some of the stars of the satellite dwarf Sph/Irr galaxy 
NGC 4627, that are shown on the separate plot 
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Fig. 2. The [{V-I), I] color-magnitude diagram of all stars located in the ACS/WFC fields of the 
galaxies. The dashed line shows the position of the TRGB and the solid lines show the regions 
of the "blue", and AGB stars. 



Holmbers j 1937 1 included 



Anj ( 1 19661) added 



show systematic dependence from extraplanar distance. The CMDs of the inner-disk areas show 
a plethora of blue and red supergiants, main-sequence stars, and blue-loop stars. The outermost 
extraplanar stellar populations are dominated by the old stars of the red giant branch. The dashed 
line indicates the position of TRGB (the tip of the red giant branch) and above TRGB some of 
the stars are asymptotic giant branch (AGB) and red supergiant (RSG) stars. 

One of the ACS/WFC fields around NGC 4631, additionally i ncludes some of st ars of the 
neighboring NGC 4627. Relying on their similar radial velocities 
NGC 4631/NGC 4627 in his catalog as a galaxy pair Holmberg 442, while 
them in his catalog of pecuhar galaxies as Arp 28 1 . The dwarf companion NGC 4627 is usually 
classified as dwarf elliptical E4,pec (NED). The increasing overdensity of the stars that lean 
towards NGC 4627 allow us to separate them from the NGC 4631 thick disk and halo, and 
examine their appearance in the CMD presented in Fig.|3l The RGB and AGB stars show a well 
noticeable feature, that complemented by some main-sequence stars at -0.4 < (V - /) < 0.3 and 
26 < I < 26, and by a few blue supergiant candidates. If even a small part of these blue young 
stars belongs to NGC 4627, the classification of this galaxy should be altered to a transitional 
dSph/Irr type. The relatively blue colour of NGC 4627 (B - V) = 0.62 (RC3), and clear evidence 
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Fig. 3. The [(V - I), I] color-magnitude diagram of the stars located above the thin disk of the 
galaxies. The dashed line shows the position of the TRGB and the solid lines show the region of 
the RGB stars. 



of the O ii-emission around its circumnuclear area JBettoni & BusonligSTh are also indicative of 
the current presence of some star formation in this galaxy. 

For our spatial structure analysis we used only the stars with magnitudes above the 50% 
completeness level (i.e. the magnitude where only half the simulated stars are identified), that 
occurs at about 27.3 < / < 27.5 above the limiting / ~ 28 mag in these diagrams. 



3.2. Distances 

The high quality of the presented data provides a valuable opportunity to estimate distances of 

investigated galax ies, based on on e the most successful and reliable distance candle, the luminos- 
I II 1 

ity of the TRGB (ILee et alJll993^ . While being affected by contamination of bright asymptotic 
giant branch (AGB) stars, statistical (Poisson) noise of star counts, and dependence from the 
metallicity/age of RGB stars, the TRGB method gives an estimate of distance with precision an d 



Bellazzini etalJ2001). 



accuracy similar to that of the Cepheids method, < 10% ( ILee et aiTlQQS 

In order to estimate distances, we have used the outer extraplanar areas in each galaxy, avoid- 
ing crowded star formation regions in a thin disk. Since fraction of AGB-to-RGB stars steeply 
declines with extraplanar distance, we can also expect of decreasing of the contamination by 
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Fig. 4. Smoothed luminosity function (solid lines) and edge-detection Sobel-filter output (dotted 
lines) for the red stars. The position of the TRGB corresponds to the peak of the Sobel-filter. 



AGB stars in the outermost areas. At the same time, we maintained the coverage large enough to 
include a statistically significant number of RGB stars. Following the standard technique we cal- 
culated the luminosity function of candidate RGB stars, and applied an edge detection algorithm 
(Sobel filter) to detect the location of the TRGB. The position of the first peak in the convolved 
function was used to determine the tip magnitude while the width of the peak provides an esti- 
mate of the tip error. The stellar luminosity functions and the results of their convolution with 
Sobel filter are presented in Fig.Bl and reveal a sudden discontinuity which corresponds to the 



TRGB. Using the method of 



Lee et al 



1 1993h . we estimate distances to the galaxies and mean 



metallicity of their red giants. The inferred metallicity of RGB stars should be consi dered solely 



as a l ocal value, due to the metallicity gradient commonly seen in spiral galaxies ( iTiede et al. 



20041) . We obtain the following results for the galaxies. 



3.2.1. IC 2233 



The TRGB in IC 2233 is discerned at /trgb = 26.09 + . 08, an d the median colour of RGB stars 



(V - /)-3.5,o is 1™60. Using the calibration bv lLee et al. 



1 19931) . we obtain a metallicity [Fe/H]^ 



-0.93 dex and a distance modulus (m - M)o - 30.09, corresponding to D = 10.42 + 0.38 Mpc. 
Note that the error in distance includes random (residuals in the PSF fitting, tip measurement. 
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extinction) and systematic (RR Lyrae distance scale, absolute magnitude of TRGB, photometric 
zero points) uncertainties. 



3.2.2. IC 5052 

Based on the luminosity function of RGB stars in the outer area of IC 5052 (Fig. 2}, the TRGB 
is at /trgb.o = 24.74 + 0.07. The median colour between /q = 25™ 1 and /q = 25™3, which we use 
to determine (V - /)-3.5,o, is 1™66. These results correspond to a metallicity [Fe/H]^ -0.9, and a 
distance (m - M)o = 28.75 (D = 5.62 + 0.20 Mpc). 



3.2.3. NGC 4631 & NGC 4627 

We have performed single-star photometry for two fields in the vicinity of NGC 463 1 . One of the 
fields (SI) partly covers the circumnuclear area of NGC 4631 and the outer part of the satellite 
dwarf Sph/Irr galaxy NGC 4627, while another area (S2) is at the larger galactocentric distance 
along the disk plane (see^J. To avoid bright AGB stars of the NGC 4631 disk and bulge, we 
selected stars of the north-east sector of the S 1 field and south part of the S2 field (marked by 
arrows). The location of the TRGB is found by detecting a sharp edge in the luminosity function 
of the red plume at /trgb.o = 25.20 + 0.05 and median colour of the RGB (V - /)-3.5,o is 1 -60 (for 
stars of SI field), while for stars of S2 field /trgb.o is at 25.17 + 0.05, and (V - /)-3.5,o = 1.47. 
These results correspond to a metallicity [Fe/H]^ -0.93 (for SI) and [Fe/H]^ -1.31 (for S2), 
and a distance (m - M)o = 29.29 (SI) and 29.22 (S2). The distance modulus of NGC 4631, 
based on the average results for two fields is (m - M)o = 29.26 + 0.04, that correspond to a 
distance D = 7.11 + 0.13 Mpc. Roughly selecting stars around the neighboring NGC 4627, we 
found the TRGB edge of their luminosity function at /trgb.o - 25.12 and the RGB colour of 
(V-/)-3.5,() = 1.38. The derived distance of NGC 4627 (m-M)o ^ 29.13 (D^ 6.70 + 0.15 Mpc) 
is less certain than the distance estimation of NGC 463 1 , due to the small star number counts and 
unclear contamination by the NGC 463 1 stars, but it confirms their proximity to each other 



3.2.4. NGC 5023 

The / - band luminosity function of RGB stars in the outer area of IC 5023 (Fig. 0}, has a 
discontinuity at /trgb.o - 24.95 + 0.05 corresponding to the TRGB edge. We used the median 
colour between /o = 25™3 and /q - 25™6 to determine {V - /)-3.5,o, is 1"35. These results 
correspond to a metallicity [Fe/H]^ -1.7, and a distance (m-M)o = 28.94(D - 6.14+0.15 Mpc). 



3.2.5. Comparison to previous distance estimations 

In the context of distance verification, it is a valuable opportunity to compare our TRGB dis 



tances estimations with a results of other methods. We also list distances calculated by 



Seth et al 



1 2005 ah as /)sdj2005, who used the similar TRGB method from HST data. 
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There are three published dista nce estimations f or IC 2233: 7.7 Mpc ( iRossa et al 



6.9 Mpc llWilcots & Prescotl 



2004 . and 10.6 Mpc iITuIIv 



D = 10.4 Mpc, does agree well with the 
from two other results. 



f 



2004 . 



1988h . Our distance measurement, 



Tullvl (Il98i^ determination, but significantly differs 



by 



For IC 5052, we found the foll owing previous estirnation s: A distance of 6.0 Mpc reporte d 



Seth et al 



1 2005al). 6.3 Mpc b v 



Karachentsev et al 



5.9 Mpc dRvan- Weber et alJ2003h . and 9.2 Mpc dBecker et al 



1 2004 



, 7.9 M pc by 



198 



Rossa et al 



1 2004 . 



While our distance for NGC 463 1 , D = 7.1 Mpc^ is close to the 



measurement of 7.5 Mpc and 7.6 Mpc by 



Seth et al 



estimation of 5.3 Mpc by 



LSofuefcHandaMl99i 



3k 



Our value is Z) = 5.6 Mpc. 



Rand & van der HulstI ( Il993 1 



1 2005 ah . it significantly disaccords with the 



There are a few distance estimations for NGC 5023. 



van der Kruit & Searld ( 119821) accepted 



a distance of 8 Mpc b ased on its re dshift and Hp = 75 km s ' Mpc Distances of 5.6 Mpc 
I Bottinelhet al.ll985 1 and 4.8 Mpc ( Swaters et ^iL ' 



il.l2002h were calculated using the Tully-Fisher 



relation. A brightest star method measured by 



Sharina et al 



( 19991) gives D = 5.4 Mpc. Ou r 



TRGB distance is 6.14 Mpc, that is in an agreement with Dsdjioos = 6.61 by 



Seth et al. 



f 



12005a I. 



4. Stellar spatial distributions 



The regions outlined in Fig. |2l and |3] indicate the regions of the AGE, "blue plum" and RGB 
stars used for the analysis of their extraplanar di stribution. Similarly to our previous studies of 



the outer stellar surroundings of th e disk galaxies dDrozdovskv et alJ2003 



Tikhonov et alJ2005 



Tikhonov & Galazutdinoval2005l) . we apply the following method. Using colour-magnitude dia- 



grams, we separate different stellar populations and analyze their spatial distribution. By sorting 
stars on the different groups according to their life expectancies, we assign probable membership 
of stars to the different galaxy components. The measurements of the stellar number density per- 
pendicular to the galaxy plane and search for the truncation radius (edge) of the thin and thick 
disks allow us to trace farther a presence of stellar halo. Unless otherwise stated, it is assumed 
that stellar thick disk and halo are close to axially symmetric, allowing us to extrapolate the spa- 
tial distribution from mapping only a part of the galaxy. Forthcoming, we will examine the case 
of two neighboring galaxies, where strong tidal forces might caused an asymmetry of the stellar 
structures. 

The extraplanar distribution of the stars (see Fig.|5}, shows a decline toward the mid-plane 
of the thin disks, due to the increase of stellar crowding and concentration of dust, rapid growth 
to some extraplanar distance, followed by steady exponential decline. The radially truncated 
exponential disks are rou tinely observed in the disk-dominated edge-on galaxies, reiving on the 



surface photometry (e.g.. 



de Griis & van der Kruiil996 ; 



and stellar number counts dTikhonov & Galazutdinovall2005l) . The surface number density pro 



de Griis et al.l2001 



Pohlen et alJ2004 ) 



files such as ones presented in Fig. |5l show, that distribution of different stellar populations has 
quite different scale-lengths and truncation radii. 
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Fig. 5. The surface number density distribution (SN) of RGB (filled squares), AGB (dots) and 
blue stars {empty squares) along the extraplanarZ-axis. 

The young blue stars reside mostly in the thin disk, that has a large concentration of dust. 
However, the majority of these stars are bright enough to be detected even near the mid-plane 
of a thin disk, and therefore they do not show the same drop-out as less luminous evolved stars 
(AGB and RGB). 



4.1. IC2233 

The intermediate-age AGB stars are generally more extended than blue stars and their extraplanar 
distribution roughly follows the exponential law. At extraplanar distance of Z ^ 1.6 Mpc their 
number density drops to almost zero. The underdensity of the AGB stars make it difficult to 
obtain representative information on their outermost extend. In our previous studies, we found 
that in galaxies with larger sample of AGB stars, such as in NGC 891, there are two gradients 
of the AGB stars distribution: the steep decline of the inner AGB components until 6 kpc, and 
the sh allow decline of the outer AGB c omponent traced at least up to 11 kpc above the mid- 
plane I Tikhonov & Galazutdinoval2005 ). The most numerous among detected, RGB stars have 
the shallowest of all stellar number density gradients. The drop of their number near the galactic 
mid-plane is due to their low luminosity. At the extraplanar distance of Z = 1.9 kpc the RGB 
stars show a break in the exponential decline of their surface density, that likely corresponds to 
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a vertical (extraplanar) edge of the thick disk. Farther above the thick disk, a density of the RGB 
stars declines with a shallower slope. By fitting the exponential law, we extrapolated stellar halo 
out to Z = 4.8 kpc. Unless there is an additional outermost stellar component not covered by the 
ACS/WFC field, the relatively small extent of the stellar halo in comparison with apparent large 
axis length (planar extent) indicates that stellar halo might be represented by oblate ellipsoid. 



Stil & Israeli ( l2002l) believed that galaxies IC 2233 and NGC 2537 are a physical pair. 
However, we did not find any apparent tidal deformations of the IC 2233 thick disk, as we de- 
tected in the interacting NGC 463 1 . Therefore, we assume an axial symmetry for the disk and 
halo components of IC 2233. 

It is valuable to compare the extraplanar size of the outermost stellar and gaseous structures 
of IC 2233. It is generally assumed that the size of the neutral hydr ogen envelopes of the disk- 
Martini 1998 : 



domin ated galaxies is several times larger than the stellar halos (e.g.. 



Swaters et al, 



20021) . Rescaled to our estimation of the IC 2233 distance, D - 10.4 Mp c, the extraplanar he ight 



Z(H l) - 2.4 kpc accordin g to measurements of 



LStil fc Israeli ll2002ll. or 



of the H I envelope is 

6.1 kpc according to fWilcots & PrescottI (120041) . We detected the edge of the IC 2233 thick 
disk at extraplanar distance of 3.8 kpc, while its stellar halo is traced as far as 9.6 kpc. Note 
that the extraplanar emis sion-line gas was traced out to ~ 1 kpc from the IC 2233 mid-plane 



I Miller &Veilleuxl2003). 



I 



4.2. IC 5052 

Similarly to IC 2233, the young blue stars of IC 5052 are concentrated in the thin disk, that 
truncated at ~ 1 kpc above the mid-plane. The intermediate-age AGB stars are confined in both 
the thin and thick disks and can be seen out to Z = 1.6 kpc. The old/intermediate-age RGB 
stars reveal break in their surface density at Z = 1.9 kpc, but extend farther with smaller surface 
density gradient tracing a stellar halo. The extrapolated vertical size of the halo is Z ^ 4.0 kpc. 



There is a reported detection of a bright H i-halo and extended emission jRossa & Dettmai 



2003 ), but the exact s ize is unknown. In the presented Hi column density contours 



Rvan- Weber et alJl2003l) . the IC 5052 vertical structure is seen at least until Z = 3.2 kpc, but 



there are also indication on more extended hydrogen filaments. 



4.3. NGC 4631/4627 

Since one of the ACS/WFC fields (SI) covers partly the circumnuclear area of NGC 4631 while 
second one (S2) are well off-center along the disk-plane, we can easily discern the bulge stars 
from the disk and halo ones (see Fig.[0. The distribution of stars detected in the field SI are 
symmetric, and show a significant extent as for the young blue stars (truncation height Z ^ 
2.0 kpc), intermediate-age AGB stars (Z ^ 4.0 kpc), and RGB stars, which spread farther than 
the SI field borders. 
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The surface density profile of stars detected in the S2 field show quite different distribu- 
tion. Only profile of blue stars is symmetric relatively to the mid-plane, while AGB and RGB 
stars show asymmetric breaks at opposite sides of galactic plane: the north side is clearly trun- 
cated at smaller extraplanar distance than the south one. For the AGB stars the South-to-North 
sides asymmetry of this break is Zsouth /Z^onhiAGB) - 2.8 4.0 kpc, and for the RGB stars 
Zsouth/ZNorth(-RGB) = 3.1- > 4.5 kpc. 

The pecuUar distribution of the evolved AGB and RGB stars can be explained by different 
possibilities: 

i) The interaction between NGC 463 1 and NGC 4627 distorted closest (Northern) side of the 
tMck disk and halo, that resulted in their spread out toward North part from galaxy plane. The 
spatial stretch of stellar structures yielded to the diminishing of their surface density. Certainly, 
some of the outermost stars in that area might actually belong to NGC 4627. Another difficulty 
of this model is to explain the high-degree symmetry of the distribution of blue stars relatively to 
the galaxy plane, as it was found in both of the studied fields, SI and S2. 

ii) The asymmetry in the distribution of evolved stars is purely due to the mixing of the NGC 4627 
stars with NGC 463 1 . This model is not supported by the data, since the surface overdensity is 
detected toward the opposite (Southern) side across the NGC 463 1 plane. 

Hi) The asymmetry in detection of the AGB stars of NGC 463 1 is due to the presence of the 
foreground dust clouds in NGC 4627, that significantly dim light of the NGC 463 1 stars, shifting 
them below the magnitude of the TRGB. This mechanism might yields to the observable un- 
dercounting of AGB stars by 2.3 times, and requires an average interstellar e xtinction of about 
one magnitude in /-b and, that is possible for ~ 1 kpc extraplanar distance iGoUa et alJ 



Martin & Kenj2001l) . The RGB stars are more spread out and thus less affected by dust extinc- 



tion. Therefore their South-to-North surface densities differ less (the ratio is 3 -H 2). This model 
is supported b y the H i-feature (spur 4), plausibly assoc iated with NGC 4627 and apparent on the 



radio-images llRand & van der Hulstll993 



tll993[|Randl 



lRandl1994D. 



We estimated the outermost vertical extension of the NGC 463 1 stella r halo as » 8 h- 10 kpc. 
Interestingly, from the spectroscopic observations iMartin & KernI ( 120011) have detected the ex- 
traplanar emission line gas around NGC 463 1 out to Z ^ 7 kpc. 



4.4. NGC 5023 

As in the case of IC 2233 and IC 5052, the stellar surface density of NGC 5023 stars are sym- 
metric relatively to a disk plane within the observed field. The young blue stars are prominent 
but confined to the narrow plane (Z = +0.8 kpc) of the galaxy thin disk. The extraplanar size 
of the disk-like structure of AGB stars is 1.3 kpc, while thick disk of the RGB stars show break 
at Z = 1.6 kpc. The halo are well traced above by RGB stars until Z - 3.2 kpc. NGC 5023 
is the most flattened among analyzed edge-on galaxies with length-to-height ratio of a/b - 9.3 
(NED). Interestingly, the stellar halo of this galaxy indicate also of the largest flatness. Does it 
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indicate the rotation of the halo? The corrugation of the NGC 5023 disk have been reported by 



Florido et al 



1 19911) . similarly to the disk of the Milky Way. The H i-disk of the NGC 5023 show 
a sharp edge at Z = 1 .8 kpc, while some of the hydrogen filaments are seen at larger distances 
S waters et al.l2002 ). 



5. Extraplanar colour gradients of the RGB stars 

Metal abundances act as a tracer to the conjoint action of star formation and dynamical mixing of 
stars and gas within a galaxy. Therefore spatial abundance variations are among the most impor- 
tant parameters in any theory of galactic chemodynamical evolution. The majority of the investi- 
gations concentrated predominantly on the radial (in-plane) abundance gradients in the galactic 
disks, which have been established in the Galaxy and other nearby disk galaxies based on the 



spectroscopy of a variety (mostly c 
nebulae, and brightest stars (e.g.. 


)f young) obiects such as stellar clusters, Hn regions, planetary 


Shaver etal.1 1983 


Friel&JaneJl991 


Rolleston et allboOO 


Marmiez etalJ2002 


Andrievskv et al.l2002l. Edse-on galaxies provide an opportunitv to explore 



the abundance gradient as a function of extraplanar height, that ensure minimal contamination by 
stars belonging to the dominant galactic disk component. However, only a few edge-on galaxies 
are close enough to perform metalhc ity measurements from the spe ctrum of evolved stellar pop- 
ulations due to their low luminosity ( Reitzel & Guhathakurt J2OO2 ). The idea to use the ( V - I) 



colour of the RGB stars comes from its much greater sensitivity to metallicity than it has to age, 
makin g the mean RGB colour a good stellar metallicity indicator, although there is some degen- 



eracy llLee et al.ll993h . The method is based on the statistical approach and its accuracy depends 
on the quality of the photometric data as well as on the number of detected RGB stars (including 
the source confusion-Galactic red dwarfs, etc.). The advantage of this method is that red giants 
can be situated at relatively large galactocentric distances, where surface brightness is very low 
and spectral methods are below the sensitivity limits.^ 

The number density of RGB stars varies several orders of magnitude across the extraplanar 
space. To include statistically significant data, we track colour changes of RGB stars located in 
three major areas: the equatorial (thin disk) zone, thick disk and halo. We consider only two 
galaxies for this study, NGC 5023 and IC 5052. NGC 4631 shows the apparent distortion of its 
disk, and make it difficult to separate stars belonging to different spatial substructures. IC 2233 
is too distant to collect enough stars for the reliable statistical study. 

Extraplanar colour distribution of RGB stars are similar in NGC 5023 and IC 5052 (see 
Fig. |6j and characterized by weak negative median (V - I) gradient out to the "halo" area. 



^ Note, that the metallicity gradients of young and evolved stars reflects different structures, and are the 
subject of the different star forming and kinematic effects: / These estimations deal with different (spatial?) 
substructures - thick and thin disks, ii The metallicity of young stars reflects rather the SFR and CEL, 
while the distribution of old stars are more depends on the effectiveness of the process of stellar mixing. We 
consider here only the metallicity of the upper (brightest) part of the RGB (i.e. intermediate-age/old stellar 
populations). 
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Fig. 6. The distribution of the (V - 1) color of the RGB stars from 0.3 to 0.7 mag below the tip of 
the RGB for three extraplanar zones: thin and thick disks (solid lines) and halo (dotted line). The 
similar morphology gradient in the RGB is visible in both of the edge-on galaxies, NGC 5023 and 
IC 5052: the red RGB stars are more concentrated to the disk plane, while outermost extraplanar 
stars are systematically bluer 



Zero or slightly negative vertical colour gradient have also been found in some other 



ously studied edge-on disk galax ies, such as NGC 55 and NGC 4244 llTikhonov et al 



Tikhonov & Galazutdinova 



previ- 



2m 



2005h . There i s no apparent changes in colour gradient within the 



thick disk area based on our measurements 



Mouldl ( l2005h found the zero or slightly positive ver- 



tical colour gradients for the RGB stars in the thick disks of four others local edge-on galaxies, 
but no data was av ailable for the z > 2 kpc area of these galaxies. Based on the similar HST data. 



Seth et al 



1 2005bl) confirmed that extraplanar stars show little or no colour gradient. 
Note, that no correction has been made for interstellar reddening, that might mostly affect the 
colours of the thin disk stars. Additionally due to the crowding, the sample of objects within a 



few hundreds parsecs of the disk mid-plane is likely dominated by the outer-disk stars with large 
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Fig. 7. A scaled 3-D represantation of the stellar components of a typical spiral galaxy and results 
of their projection as a stellar number density for a face-on (bottom plot) and edge-on (right) 
galaxy. 



galactocentric radii. The lack of the inner-disk RGB stars, as well as age-metallicity bias might 
significantly impact the cumulative colour distribution. Thus the lack of strong colour gradient 
do not ruled out the systematic chemical changes of RGB stars perpendicular to a disk plane, 
as well as star to star variations in elemental abundance ratios. Even the outermost "halo" RGB 
stars apparently show a large range of colours, indicating that extraplanar areas have experienced 
a complicated star formation history with various self-enrichment processes. 



6. Model of the stellar structures of a disk galaxy 



Comparing the inferred parameters for the thin & thick disk and halo with our and other published 
results, we compose a scale 3-D model of the stellar structures of a typical disk galaxy (Fig. 01. 
This model relies on the star number counts in the low-inclination, M 8 1 and NGC 300, and 
eight high-inclination disk galax ies, NGC 55, NGC 891, NGC 4144, NGC 4244, NGC 4 363, 



NGC 5023, IC 2233, IC 5052 llTikhonov et al 



Tikhonov & Galazutdinova 



2m 



this 



paper). Now we have new results of stellar photometry of four disk galaxies: NGC 1311, 
UGC 1281, UGC 8760 and IC 1959 that to support this model of the stellar structure. While 
our primary goal in this work is to constrain the spatial geom etry of the spiral galaxies , we would 



like to note that massive irregular galaxies, such as IC 10 (iDrozdovskv et al 



20031) and M 82 



Tikhonovl2005ft . show all kinds of the stellar structures, including thick disk and halo. The spa- 
tial distribution of stars along and perpendicular to a disk galaxy plane are alike. The young blue 
stars are mostly confined to the thin disk of the galaxy, while elder stars are more spread out with 
oldest stars are traced outermost. 
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Table 2. Properties of Galaxies 



Name R.A.(2000.0) DEC(2000.0) 



axb 



B, 



Type 



A, 



m- M 



Mb 



1C2233 OSHS^SS'igi +45°44'31'.'7 565 5.2x0.6 13.07 11.70 SB(s)d:sp 0.171 90 30.09 

IC5052 20''52"'0H63 -69°11'35'.'9 598 5.9x0.8 11.16 10.52 SBd:sp 0.168 90 28.75 

NGC4631 12''42"'08'!01 +32°32'29'.'4 606 15.5x2.7 9.75 8.61 SB(s)d 0.056 90 29.26 

NGC5023 13''12"'12':60 +44°02'28'.'4 407 7.3x0.8 12.85 11.59 Scd 0.060 90 28.94 



-18.39 
-18.23 
-20.65 
-17.35 



The Galactic extinction correction is by 
The inclination is taken from LEDA. 
The value of (m - M) is from this paper. 



Schlegel et al. 



199^) 



7. Discussion 

Relying on the star number counts in the disk galaxies IC 2233, IC 5052, NGC 4631 and 
NGC 5023, we examine the extraplanar spatial distribution of young and evolved stars in their 
disks and halos. These edge-on disk galaxies show the similar morphological properties, thin and 
thick disks and halo, but relative sizes of these structures are different in each galaxy. Combinin g 



this work with re sults obtained for other four large edge-on disk galaxies 



Tikhonov et al. 



1 20051) : 



Tikhonovl ll2005l) . we suggest that most of the massive disk galaxies have not only a thick disk 
but also an extended halo, consisting mainly of evolved stars. 

The stellar number density of the halo at the truncation radius (break) of the thick disk varies 
a lot, and plausibly correlate with overall flatness of a disk galaxy. The ratio of the halo-to-thick 
disk size is relatively constant, ^haio/^thickdisk - 2.5 + 0.5. It might reflect the overall rotation of 
the galaxy, including flatten stellar halo. 

Noteworthy the significant contrast between stellar structures of two nearest giant spiral 
galaxies, the Milky Way and M 31, the only systems there possible to measure the kinemat- 
ics and chemical properties of individual stars. While both of the "sister" galaxies are about 
the same baryonic-plus-dark masses, their stellar disks and halo differ significantly. The outer 
second disk-like stellar structure of M 31 are 3-4 times more extended than in our Galaxy and 
domi nate over the halo stars at least out to a de-projected galactocentric radii of about 40 kpc 
Worfliev et al.ll2005t 



(e.g.. 



Ibata et al.ll2005l) . The significant difference in the properties of the 
outer stellar populations of the nearby disk galaxies, such as the Milky Way, M 31, and others 
(including our dataset) plausibly suggest that the final morphology of a galaxy is extremely sen- 
sitive to the details of the satellite accretion history, as well as peculiarities of the star formation 
and feedback. 

Such unique populations as bulge, thin and thick disks, and halo have likely formed by differ- 
ent mechanisms, especially considering the differences in ages and metallicities. The thick disk 
and halo may cast light on the most ancient structure formation, so being extremely important 
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for the understanding of early history of galaxy formation and evolution. If stellar halo have been 
formed prior to the thick disk (see a recent review of various formation scenarios for the stellar 

should show both the characteristics age and metallicity pr ofile of 



halo by 



a pos thalo epoch of star formation 6-12 Gyr ago and differ significantly from the halo ( iMoulc 



20051) . The addition of results from photometry of more evolved stars (e.g., Horizontal Branch 



and old Main-Sequence Turn-Off stars), as well as from kinematic and chemistry of resolved 
stars are required to refine various formation scenarios. 



8. Summary 

On the basis of the star count method, we have analyzed extraplanar stellar structures in four 
large high-inclination galaxies. The main results can be summarized as follows: 

- The extended stellar thick disks and halos have been detected in IC 2233, IC 5052, NGC 463 1 
NGC 5023. 

- There are clear differences between surface density gradients of evolved stellar populations 
assigned to the thick disk and halo of these disk galaxies, which allowed us to detect the edge 
of the thick disk. The extraplanar distribution of the outermost stars indicate that stellar halo 
is flattened along the extraplanar direction. 

- Having a large, high quality photometry of the RGB stars allowed us to estimated the dis- 
tances to the galaxies from the TRGB magnitude with high confidence. We determine a 
distance of 10.42 + 0.38 Mpc for IC 2233, 5.62 + 0.20 Mpc for IC 5052, 7.11 + 0.13 Mpc for 
NGC 4631 (6.70 + 0.15 for its satellite NGC 4627) and 6.14 + 0.15 Mpc for NGC 5023. 

- In NGC 4631, the revealed asymmetry of its stellar thick disk and halo, plausibly is an evi- 
dence of an past ore ongoing interaction with NGC 4627. 

- The stellar distribution of stellar populations of these four galaxies provide more confidence 
for the unified model of three-dimensional structure of large disk galaxies, and motivate 
follow-up spectroscopic study to constrain the kinematic characteristics of the stellar compo- 
nents. 
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